THE IDENTIFICATION PROBLEM FOR THE ATTENUATED X-RAY 

TRANSFORM 



PLAMEN STEFANOV 

Abstract. We study the problem of recovery both the attenuation a and the source / in the 
attenuated X-ray transform in the plane. We study the linearization as well. It turns out that 
there are natural Hamiltonian flow that determines which singularities we can recover. If the 
perturbations 5a, 8f are supported in a compact set that is non-trapping for that flow, then the 
problem is well posed. Otherwise, it may not be, and least in the case of radial a, /, it is not. We 
present uniqueness and non-uniqueness results for both the linearized and the non-linear problem; 
as well as a Holder stability estimate. 



1. Introduction 
We study the attenuated X-ray transform 

(1.1) X a f(x, 0) = j e -Ba(x+te,e)f( x + ^ dt) x e R 2 , # e S" 1 , 

in the plane with a source / and an attenuation a that we want to recover. We denote by 

(1.2) Ba(x,6) = / a(x + t9)dt 

Jo 

the "beam transform" of a, usually denoted by Da. We will assume that both a and / are compactly 
supported. In applications, a constant attenuation a is also considered but when observations are 
made on the boundary of a compact domain, one can replace that constant by a constant multiple 
of the characteristic function of that domain. 

The problem that we study is: can we recover both a and / from knowledge of X a fl Sometimes 
this is called the Identification Problem (for SPECT). 

This problem arises in Single Photon Emission Computerized Tomography (SPECT). Radioactive 
markers are injected into a patient's body and the emitted X-rays, attenuated by the body, are 
detected outside of it. The problem is to recover the source with a unknown attenuated coefficient. 

When a is known, it is well known that / can be reconstructed uniquely, even by means of explicit 
formulas (6j[T], [221 123| . |21j . For this reason, some of the numerical attempts to do a reconstruction 
are focused on recovery, or getting a good approximation of a first, instead of treating (a, /) as 
a pair. Sometimes this is called attenuation correction, see e.g., [MIEl]. In clinical applications, 
additional X-rays are taken to reconstruct a first. Eliminating or reducing those additional X-rays 
remains an important problem. 

There has not been much progress in the mathematical understanding of the identification prob- 
lem so far. A related but not identical problem for finding both a constant attenuation and the 
source in the exponential X-ray transform has been solved in |25| , see also [13] . The main result in 
[23] is, roughly speaking, that specific pairs of constant a and radial / cannot be distinguished but 
all other pairs can. The identification problem with / a finite sum of delta sources has been studied 
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in [181 H5] , see also [3], but the results there do not and cannot imply uniqueness. Natterer also 
viewed the problem as a range characterization problem: if the ranges of X ai and X a2 happen to be 
the same, for example, then there cannot be uniqueness. Range conditions, see e.g., [23], have been 
viewed as a possible tool for solving the problem, both numerically, see e.g., [5] and analytically, 
as in the recent work [2]. Numerical reconstructions have been tried, too, with variable success, in 
[H [TTl [Ml [Ml HI EH], for example. Some of them use clinical data. A. L. Bukhgeim [7] recently 
outlined a recovery algorithm if a is a priori known to be a constant multiple of the characteristic 
function of a star-shaped domain. 

Our approach is based on the following. The attenuated X-ray transform, and its linearization, 
carry information about / and a along each line twice because we integrate both in 9 and —9 
directions. From a microlocal point of view, those two lines determine the wave front sets at 
covectors normal to them. So we have two equations for two unknowns. We study first a linear 
problem that appears as a linearization of X a f near some fixed (ao,/o)- Also, the non-linear map 



A^a 2 /2 — X ai f\ is of that form, see (5.3). This problem can be formulated as the inversion of 



Zg := I Wl gi + I W2 g2, 9 = (91,92), where I w is the weighted X-ray transform with a weight w(x, 



see (2.4). The weights 1012 are of specific type in the case of the Identification Problem but we 
study general weights first. The operator X is a Fourier Integral Operator but we do not study it 
directly. Instead, to analyze the equation Ig = h, we apply an explicit operator Q to convert the 
equations 

(I Wl 9i + I W2 92) (z, ±9) = h(z, ±9) 
to equivalent pseudo-differential ones of the type 

(1.3) (wi(x, ±D ± /\D\) + l.o.t.) gi + (wi(x, ±D ± /\D\) + l.o.t.) 9l = Qh, 



see Proposition 6.1 Here, "l.o.t." stands for "lower order terms", and Wj(x, ±D /\D\) are pseudo- 
differential operators (^DOs) with symbols Wj(x, ^z^/l^l). We view this as a 2 x 2 system of ^DO 
equations. The determinant of the principal symbol of the is given by po(x,£) = W(x, C" 1 /^!), 
where 

W(x, 9) = w\{x, 9)w2(x, —9) — wi(x, —9)w2(x, 9). 

Since W is an odd function of £, po is not elliptic over any x, and has a non-trivial characteristic 
variety regardless of what W\ 2 are ! m the cotangent bundle of any domain. Then po is a Hamiltonian 
of fundamental importance for this system. The singularities of g that may not be recoverable lie 
on zero bicharacteristics of that Hamiltonian; moreover each zero bicharacteristic either consists 
of singularities only, or there is none on it. This brings us to the following condition, well known 
in the theory of ^DOs of real principal type: if g has a support in a compact set K that is non- 
trapping for the Hamiltonian flow, then the singularities of g can be recovered, with a loss of one 
derivative. Otherwise they may not be but the non-trapping condition is known not to be "if and 
only if . In sections [4] and [5] based on the microlocal understanding on the problem explained 
above, we prove actual injectivity and stability of X for / supported in non-trapping K for generic 
(w\,W2), including ones satisfying some analyticity assumptions; or for small K. We then apply 
this analysis to the non-linear Identification Problem to get local uniqueness and Holder stability in 
a neighborhood of generic (a, f) under the a priori assumption that the perturbations are supported 
in a non-trapping set. 



The microlocal consequences of (1.3) are analyzed in more detail in Section 6 In particular, we 



describe the "null eigenspace" at the characteristic points. In non-degenerate cases, (1.3) is of rank 
one on the characteristic variety po = 0. 
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We also study the case of radial a and / in Section [7j A thorough study of the radial case is 
behind the scope of this work however. The reason we include it is to present an example where the 
Hamiltonian flow can be explicitly computed. The projections of the zero bicharacteristics happen 
to be the circles \x\ = R, R > 0. Then K C R 2 is non-trapping if and only if it does not contain an 
entire circles centered at 0, including the origin. In case K is trapping, and contains a ball |x| < R, 
then the uniqueness fails and there is an infinite dimensional family of pairs (a, /) with the same 
data. They consist of radial a and /. This fact agrees with the microlocal analysis that we present 
because the latter implies we may not be able to recover radial singularities. In this case actually, 
the non-trapping condition is also necessary for the problem to be well posed. 



2. Preliminaries 

The attenuated X-ray transform results from the following transport equation model. Let f(x) 
be a compactly supported source of particles (or a signal propagating along lines with unit speed) 
propagating in a medium with attenuation coefficient a(x). Then at the point x G R n and direction 
9 £ S*" -1 (the dimension n can be arbitrary), the total number of u(x,6) of particles originating 
from the source solves the transport equation 

(2.1) (9 ■ d x + a)u = f, u\e. x ^ = 0. 

This is a linear ODE along the lines t \— > (x + t9, 6) and its solution is given by 



(2.2) u(x,9) = f e-f?< x+s V ds f(x + tO)dt. 

J — oo 



This formula can be interpreted as the superposition of all attenuated signals at (x, 9) coming from 
the source. Then at points x so that 9 ■ x 3> 0, one has u = X a f. 

It is useful to extend the definition of B, see (1.2), to functions / depending on both x and 9: 

/•oo 

(2.3) Bf(x,9)= f(x + te,9)dt. 



For such /, the solution to ( 2.1 ) is given by ( 2.2 ) again, with f(x + 19) replaced by f(x + t9, 9). 
We introduce also the notation 

(2.4) I w f(x,9) = J w(x + t9,9)f(x + t9)dt, 

for the weighted X-ray transform with weig ht w{x,9). Then I w = X a for w = e~ Ba but we will 
allow more general weights in Iw 

Also, h = Xq. 

We will also denote 

v L := (— 1>2, vi), j) = (i)i,t) 2 )eR 2 . 

2.1. A Radon transform type of parameterization of X a and Since for a fixed direction 
9, x and x + s9 parametrize the same (directed) line, we will think of X a f and I w f as parameterized 
by (z, 9), z € We denote by Z the variety 

Z = {(z,9); 6eS 1 ,ze 6 L }, 

which is essentially the tangent bundle of S 1 . Then we can set z = p9 ± , and write X a f as 

(2.5) X a f( P 9 ± ,9) = J e -Ba( P e^+te,e) f ( pe ± + ^ dt> ^ 9)£¥LxS\ 
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and similarly for I w f. We think of (p, 9) £ R x S 1 as a parameterization of Z. We also define a 
measure on Z by dz := dpd9, where d9 is the natural measure on S 1 given by d$, with i? being 
the polar angle of 9. 

2.2. Functional spaces. We will assume throughout the paper that supp/ is contained in a fixed 
compact set; and we can always assume that this compact set is included in (— ir, it) 2 . We can 
therefore assume that / is defined on the torus T 2 . For any compact set if C T 2 , we define H S (K) 
to be the closed subspace of H S (T 2 ) of functions supported in if. In other words, the Sobolev 
norm in K is defined through Fourier series. We define the Sobolev spaces H S {Z) in a similar way. 



Since |p| < it in (2.5), we can assume that p belongs to the unit circle represented by [— vr,7r] with 
both ends identified. Then (p, 9) 6 Sp x Sq. The space H S (Z) is then defined by the norm 

(2-6) \\9\\Hs(z) = \\(l-%) s/2 g\\ L 2 {zy 

where d 2 is the second derivative w.r.t. p on the compact manifold S . Notice that there are no 
9 derivatives in this definition, see also [2Q|, Theorem II. 5. 2] for involving the 9 derivatives when 
a = 0. In other words, H S (Z) is defined through Fourier Series in the p variable. 

3. Linearization 

We are going to compute the linearization of the identification problem starting from formula 



(1.1). Another way to do this, based on the transport equation, is presented in section |5j 

Assume that a and / are smooth enough so that the calculations below make sense. Denote by 
G = 9 ■ d x be the generator of the geodesic flow on TR 2 w.r.t. the Euclidean metric. Since a has 
compact support, then GBa = —a, and Ba = for x ■ 9 ^> 0; and Ba = I\a for x ■ 9 0. Here, 
I\ = I a for a = 1. 

Since the problem is linear w.r.t. /, we linearize near some a first, with / fixed. Let a s = a + s5a. 
Then 

A x a J = - [ e -Ba{x+te,e)f, + td \ B5a r x + w d \ dt 
ds s=o J 

Write 

e -Baj = _Q Be -Baj 

and plug this into the formula above. Integrate by parts to get 
d_ 

ds s=o 

The linearization of X a f w.r.t. a is therefore a weighted X-ray transform of the perturbation 5a of 
the form 



X a J = J [{Be- Ba f) 5a] (x + tO, 9) dt - X a f.X 5a. 



w(x + t9,9)5a(x + t9)dt 
with a weight function 

(3.1) w = Be- Ba f-X a f. 

The second term on the right is constant along each line. The weight can also be expressed as 

(3.2) w(x, 9) = - e -Ba(x+t0,e) f ( x + 19 ) dt 

J — oo 

A direct calculation yields 

(3.3) w = -e~ Ba u, 



where u is the solution (2.2) of (2.1). 
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Let 8X a j(Sa,Sf) denote the linearization of X a f near fixed a, /. We just proved the following. 
Proposition 3.1. 

5X aJ (5a,5f) = I w 5a + X a 5f, 



where w is as in (3.2) or (3.3). 



4. A MORE GENERAL LINEAR PROBLEM: INVERTIBILITY OF A SUM OF TWO WEIGHTED X-RAY 

TRANSFORMS 

4.1. Formulation and preliminaries. Consider a more general problem. Let I(gi,g 2 ) = I Wl 9i + 
Iw 2 92i where w\^ 2 are two weight functions, i.e., 

(4.1) I( gi ,g 2 )(x,9) = J w 1 (x + t9,9)g 1 (x + t9)dt + J w 2 (x + t9,9)g 2 (x + tff) dt. 

We will compute X*X w.r.t. the inner product in L 2 (Z). Clearly, 

( T* T T* T \ 

(4.2) X*X = [ t Tt r*r 

By Proposition A. 2 1*1 is a ^DO of order —1 with principal symbol 

2 



(4.3) 
where 



5a v {I*X) 



7T 



|wi,+r + wi j+ w 2l + + wi-w 2 - 

\w 2 +\ 2 + \w 2 J 2 



w jt ± = w j (x,±Z ± /\Z\), 3 = 1,2 



A direct calculation yields 

(4.4) det 5a p {X*X) = \w\^j r w 2 - — u;2,+ii'i,-| 2 



det [^+ W2 >+ 

Wl _ W 2 - 



That determinant not being zero is a microlocal ellipticity condition. As we see below, it vanishes 
over any point x; therefore, X*X cannot be elliptic over (i.e., in the cotangent bundle of) any domain. 
Set 



(4.5) 



W(x, 6) = w\(x, 9)w 2 (x, —9) — wi(x, —9)w 2 (x, 9). 



Then det Sa p (X*X) = \W(x, ^ _L )| 2 . The function W is odd in 9, and therefore, for any x it has 
zeros for some vectors 9. The inconvenience of working with (4.4) however is that it has double 
characteristics. 

Instead of studying the invertibility of X*X, we will approach the problem in a more direct way, 



slightly different (but equivalent) than what we do in Section [6j see also (1.3). Set 

(4.6) Jh(x,£) = h(x,-£). 

Let a(x, 9) be any smooth function, odd on S 1 w.r.t. 9. Apply I' ajw2 to the equation 

(4.7) I Wl gi + I W2 g 2 = h 
to get 
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By Proposition A.2[ both operators on the left are ^DOs of order —1. The principal symbol of 



I 'aJw 2 I wi is S iven b y 2n /\£\ times the even P art of (awiJw 2 )(x,£ "VICD) Le -> b y 27ra(x,^ "Vl^D/ICl 
times the odd part of (w\ Jw2)(x, ^ Thus 



(4-9) a p (I aJw J Wl ) = —-aW\ 0= £ ±/ ^. 



Notice that W is the determinant in the r.h.s. of (4.4) but not squared. It has the same zeros as 
(4.4) but they are simple. In the same way, we get that the principal symbol of I' a j W2 Iw2 is as 



above but with ui2 replaced by wi, i.e., it is zero; and therefore, l' a j W2 Iw 2 is of order —2. 

Choose a = 0\ first. Then |^|a(^ _L /|^|) = —£2, which is the symbol of —D2 = 182, and we get 

(4.10) M4j™ 2 V) = T^(-& w ( x > e )\e=^m- 



Modulo lower order terms, (4.8) becomes 

^_(-D 2 )W (x, D ± /\D\) 91 - I' eiJw2 h, 

where the meaning of 1/|D| is given by the ^DO calculus. Similarly, taking a = 82, we get 

^_D 1 w(x,D ± /\D\)g 1 ^l' d2jw2 h. 

Apply —D2 to the first identity, D\ to the second, and add them together to get 

irW(x, D ± /\D\) 91 - {-D 2 l' 01 j W2 + D x l' e2j ^) h. 

Notice that the lower order terms on the left involve 52 as well. In a similar way we get 

nW(x,D ± /\D\) g 2 = {D 2 I' eiJwi - D.I'^) h. 

We therefore proved the following. 

Proposition 4.1. For all compactly supported g\, 92 we have 

(4-11) Pg = (^rHd^j^ - d 2 I , eiJw2 ,-d 1 I , e2 j m + d2l' eiJwi )lg, 

where g = (<7i, 52); and P is a matrix valued classical ^DO of order with a scalar principal symbol 
given by 

(4.12) Po(x,0-=w[x,^m). 



We notice that (4.11) can also be written in the form 
mPg 

' (6 ± -d x )w 2 (x,-6)lg(x-{6-x)6,6)d6, [ (6 L - d x )w 1 {x,-e)Xg{x - (0 ■ x)6,6) d6 



Let 

(4.13) E = {(x, £) € T*R n \ 0; p Q (x, £) = 0} = {W = 0}^ 

be the characteristic variety of po, where the sign _L applies to the second variable 9 only. 

There are several definitions of real principal type ^DOs in the literature, including or not the 
differential condition below, or the non-trapping one, in a fixed domain. We will use the following 
one. We say that the \I/DO P £ \E' m is of real principal type, if its principal symbol p m is real, scalar, 
homogeneous in £, and dp m is not collinear to £dx on {p m = 0} for £ 7^ 0. The latter condition 
says that if we identify covectors of different length by their direction, then the Hamiltonian vector 
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field never vanishes, and in particular, the flow does not have stationary points. Such operators are 
microlocally equivalent to d Xl modulo lower order terms. We also note that this condition makes 
{W = 0} a codimension one (dimension 2) smooth submanifold. The same applies to E, considered 
as part of the unit cotangent bundle. 

Below, dg± is the angular derivative in the 9 variable, i.e., the derivative d/dfi w.r.t. to the polar 
angle 

Proposition 4.2. The ^>DO P is of real principal type in some domain £1 C R 2 , if an only if 

(4.14) W , 9 ■ d x W , d g ±W cannot be all zero at the same time, for any point in £1 x S 1 . 

Proof. Extend W to 9 ^ as a homogeneous function of order 0. We have 

dp = d x W{x, ^) dx + (d/d£)W(z, Z 1 ) d£ 

= W x (x,( ± )dx + (We 2 ,-W ei )(x^ ± )d^ 

Then P is of real principal type if an only if W(x,^ ± ) = and (d/d£)W(x, = imply that 
d x W(x, £-*") is not collinear with £. The latter is equivalent to the requirement that d x W(x,^ ± ) is 
not normal to i.e., • d x W(x, ^ 0. Set 9 = Then the requirement is the following: if 
W(x,9) = and (d e ±W)(x,9) = 0, then 9 ■ d x W(x,9) / 0. Notice next that the radial derivative 
d/d\9\ of W vanishes on W = 0, therefore the requirement on dgW is actually a requirement on 
the angular derivative only. □ 

The following condition plays a critical role in the theory of local solvability of ^DOs of real 
principal type. Let K be a compact subset of R 2 . 

Definition 4.1. We say that K is non-trapping (for po) if there is no maximally extended zero 
bicharacteristic that lies entirely over K. We call the projections of the zero bicharacteristics of po 
to the x space rays. 

Notice that the rays are continuous but not necessarily smooth curves. They may even degenerate 
to a point, see Example |4.1[ 

Hormander's propagation of singularities theorem (see, e.g., |31l VI. 2.1] or |16l Theorem 26.1.4]) 
implies that if g is supported in a non-trapping K, and Pg G H s , then g E fP -1 . In other words, 
we have non-local hypoellipticity, with a loss of one derivative. As a consequence, by the open 
mapping theorem, for any s and £, there is C > so that 

(4-15) \\g\\ HS -i < C\\Pg\\ H . + C\\g\\ H i, V 5 G C?(K), 

see eqn. (VI. 3. 3) in [51] . 

Example 4.1. Let w\ = \9 ■ x and W2 = 1. Notice that w\ is not non- vanishing. In this case, 
W = 9 ■ x. Then \£\po = x ■ = — a?i£2 and E consists of (x,£), £ 7^ that are collinear. 
In other words, all singularities that may not be recoverable are the radial ones. The Hamiltonian 
equations are given by 

Xl=X 2 , X 2 = -X\, £l=6, &2 = -£l- 

The zero bicharacteristics then are given by 

(4.16) x = i?(sint,cost), £ = A(sint, cost), R>Q,\^Q. 

Their projections on the base (the rays) are given by the circles x = R(s'mt, cost), R > 0. If R = 0, 
then that projection is a point. The whole bicharacteristic is not stationary however and is given 
by x = 0, £ = A(sint, cos t), A 7^ 0. We then see that a compact set K is non-trapping if and only 
if K contains no entire circle \x\ = R, R> (including the origin), see Figure ^1 Then Xg recovers 
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the singularities of g = (51,52)- If K is trapping, the singularities that may not be possible to be 
recovered are the radial ones. 

Inverting Xg = I Wl gi+Iig2 is easy. The first term is odd w.r.t. 9, and the second one is even. The 
equation Ig = h then decouples into two equations I Wl g\ = /i dd, I\92 = ^evcn- The kernel of X (on 
£'(R 2 )) then consists of pairs (51,0), where g\ G Ker/ Wl . Using arguments similar to the Fourier 
Slice theorem, see Section[6j we can see easily that I Wl 9\ = if and only if {x 2 d\ — x\d 2 )gi = (the 
operator in the parentheses is po(x,D), up to an elliptic factor) and the solutions of the latter in 
£'(R 2 ) are given by all compactly supported radial distributions. In particular, on L 2 omp (i?(0, 1)), 
the kernel of I m consists of all radial functions in that space. We therefore get 

(4.17) KerX = {(g u 0); g x is radial} . 

Since radial functions can have (radial) singularities at all points, we get that no radial singularity 
of <7i, i.e., a singularity of the type (x = £), /i £ R, ( / 0, can be recovered in general. On 
the other hand, if K is non-trapping for po, and supp<? C K, then they can. In that case, I Wl is 
microlocally equivalent to the derivative w.r.t. to the polar angle in R 2 ; and since on each circle 
there is an open arc where g = 0, that circle cannot support a singularity of g, by the propagation 



of the singularities theorem. In fact, by (4.17), if K is non-trapping, we have more: g± = 52 = 0. 
This example also reveals that Ig G C°° is not microlocally equivalent to (4.11), see Section [6] 



for more details. Indeed, only the radial singularities of g± are not recoverable, while those of 32 are 
recoverable. Moreover, assume that Xg G H S (Z). Then g 2 G H s - 1 / 2 by the usual inversion results. 
We can think of I wi gi as the Doppler X-ray transform of the vector field g\(x)x. It is well know 
that we can only reconstruct the curl of g±(x)x that is (x±d 2 — x 2 d\)gi, an d the latter is in H s ~ 3 ^ 2 . 
Let suppgi G K, with K non-trapping, for example, assume that the ray {xi > 0, x 2 = 0} does 
not intersect K. Then, in polar coordinates (r, </>), 



gi(r,4>)= [(xid 2 - x 2 di)g] (r, ip) dtp. 
Jo 

This integration is not smoothing (not in all directions), and we still have g\ G H s ~^/ 2 , with an 



improved regularity in angular directions. This is consistent with (4.20) below but as we see, the 
one derivative loss is only in gi, and the that estimate does not reveal the extra regularity in 
characteristic directions. The latter is however reflected by the fact that WF(gi) can have radial 
directions only. 

It is interesting to know when the rays are smooth curves. The projection of a bicharacteristic 
to the x variables, with its parameterization determined by the Hamiltonian equation, at some 
x, has a tangent vector d^po evaluated at some (x,£) G S (i.e., Po(x,(,) = 0). This projection is 
non-degenerate, and therefore, that ray is a smooth curve, if d^po 7^ 0. There might be more than 
one £ with that property but there is at least one £ (and the whole line spanned by it) because p$ 
is odd in £. On the other hand, ^ ■ d^po = on E, therefore a tangent vector is actually £ and the 
whole line that it spans. 

Translating this in terms of W, see ( |4. 12 ), we get the following. If for some (xq, #0), we have 



(4.18) W(x , 9 ) = and d 6 ±W(x , 9 Q ) ± 0, 

then there is a smooth ray through it. Moreover, starting from (xq, 9q) with that property, by the 
implicit function theorem, we can solve W(x, 9) = locally for 9. This gives us a smooth unit 
vector field, with integral curves that are rays. Then W = => d 9 ±W 7^ on K x S 1 is a sufficient 
condition for all rays through K to be smooth. 
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4.2. Basic Properties of X. Below, S'(K) stands for the space of distributions supported in K, 
and we similarly define Hq(K). Also, since g is vector valued, L 2 (K), H S (K) below are spaces of 
vector valued functions. 



Theorem 4.1. Letl be as in (4-1) with w\ andw2 smooth. Let K C R 2 be a non-trapping compact 
set. Then for any s > we have the following. 

(a) For any I, there exist a constant C > so that 

(4-19) \\9\\h-(K) < C\\lg\\ H s +3/ 2 {z) + C|| 5 ||^ (R2) , Vg G C^(K). 

(b) The kernel of I on £'{K) is finite dimensional, and consist of C^°(K) functions g. 

(c) On the orthogonal complement o/KerX in Hq(K), we have 

(4-20) \\9\\h«(k) < C\\lg\\ HS+3/ 2 {z) . 

In particular, if I is infective on C^°(K), then (4.2(fy holds. 

Proof. To prove (a), we will apply ( 4.15[ ). To this end, replace P by xPx-, where \ G Co° ec L ua l s 1 
near K (this makes P properly supported, in particular). Let VL be a bounded domain containing 
suppx. Use estimate (4.15) combined with (4.11) to get for any fixed £, 

2 

\\9\\hs(K) < c Yj W^Jw^aWH^in) + C\\g\\ H t 

< C\\lg\\ Ha+3/2{z) + C\\g\\ H e, V 5 G C%°(K), 



see also Proposition A. 3. Notice the one derivative loss in this estimate since X is of order —1/2, 
see the Appendix. If we replace X by a single weighted X-ray transform I w with a non-vanishing 
weight w, then one has the same estimate but with H l l 2 (Z). We also note that Tg has compact 
support. 

Consider (b). Every g G £'{K) in the kernel of X must be smooth by propagation of singularities 
and by the assumption that K is non-trapping. Apply then (4.19) to get 

\\g\\m(K) < C\\g\\ L 2 {K) , Vg G KerXn£'(K) = KerXn C^(K). 

Since the inclusion Hq(K) <^-» L 2 (K) is compact, we get the finite dimensionality of KerX on K. 
Consider (c). Let V be the closure of C^°(K) under the graph norm H^Hffs^ + \\Ig\\ HS+S / 2 ^ z y 

We consider now X as an operator from T> to H S+3 ^ 2 (Z). Then X is a well defined bounded operator. 
Indeed, V is a subspace of the space of the compactly distributions, together with the topology. 
Then X can be considered as an operator originally defined as X : £'(R 2 ) — > £'(Z), and then 
restricted to T>. We then get 

\\g\\v < C\\lg\\ H s+3/2 {z) + C\\g\\ H i, Vg G V. 

By (a), X is injective onPfl (KerX) -1 -. Then by |311 Proposition V.3.1], for £ < s, we have the same 
inequality as above on T> n (KerX)- 1 - but without the last term. We refer also to [2SJ Lemma 3] as 
well for similar arguments, or to inequality (26.1.6) in [16]. □ 

4.3. Conditions for injectivity of X. 

Corollary 4.1. Let w\ and W2 be smooth. Let xo G R 2 be such that 
(4.22) W(x ,9) = =>- d e ±W{x ,6) / 0, V# G S 1 . 

Then ifO < e < 1, X is injective on distributions supported in the ball B(xo,e), and in particular, 
(JJfy holds for K = B(x ,e). 
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Proof. Condition (4.22) guarantees that any ray through xq is smooth at xq, and there are finite 
number of such rays. There is Eq > so that B(xo,Eq) is non-trapping. 

Assume the opposite. Then for any e = 1/j, j > 1, there is a non-trivial Cq° function supported 
in B(xo, 1/j) in the kernel of I. Then we get an infinite number of non-trivial functions <f)j in the 
finitely dimensional space V = KerXnCo°(B(xo, £o))> see Theorem |4.1[ b), with supports shrinking 
to the point xq. This is a contradiction. Indeed, —A : V — > —AV must be a bounded operator. 
On the other hand, —A is bounded below on Hq(B{xq, 1/j)) H H 2 by its first eigenvalue fij, that 
tends to infinity as j — > oo. Therefore, (—A<j)j,(j)j)/\\(j)j\\ 2 — > oo, which is a contradiction. □ 

Theorem 4.2. Let w\, ui2 be analytic inftxS 1 , where ft is an open set containing a non-trapping 
compact set K C R 2 . Then the operator X, restricted to £'{K), is infective. 

Proof. We use a result about propagation of analytic singularities, see [IT], for analytic ^DOs 
with real principal symbols. The result in [II] covers in fact a more general class of operators with 
complex- valued principal symbols that have real bicharacteristics and carries over to operators with 
matrix lower order terms. 

The operator A is an analytic ^DO in ft of order 0. Indeed, to prove that, it is enough to 



prove that operators of the kind I' a Ib, see (4.15), are analytic ^DOs ([33hof order —1 when a, b 



are analytic in ft x S . The amplitude of such an operator is given by (A.2), and it is clearly an 
analytic one , see also the proof of \29\ Proposition 1]. 

The propagation of singularities result in [11] then implies that each zero bicharacteristic of P 
in K either consists of (analytic) singular points only, or does not intersect the analytic wave front 
set of g. Since K is non-trapping, we have the latter alternative. Therefore, the analytic wave front 
set of g is empty. Then g is analytic. Since g is of compact support, we get g = 0. □ 

4.4. Generic injectivity of I. Let K be a non-trapping compact set. Then any small enough 
compact neighborhood K' of K is still non-trapping, see the proof of [16, Theorem 26.1.7]. There- 
fore, there exists an open ft D K so that every compact subset of ft is non-trapping. Then P is of 
real principal type in ft, by the definition in |16| . 

Definition 4.2. The set ft is said to be pseudo-convex w.r.t. P , if any compact subset is non- 
trapping, and for any compact set K\ C £1, there exists a compact set K2 C ft, so that every 
bicharacteristic interval in ft having endpoints over Ki, lies entirely over K2. 

In particular, if K\ is convex w.r.t. the bicharacteritics (i.e., one can choose K2 = K\), then K\ 
is pseudo-convex. 

Pseudo-convexity is a condition that guarantees existence if a global parametrix of P, see [9] 
and |16l Theorem 26.1.14]. Under that condition, we show below that injectivity of I is preserved 
under small perturbations of the weights. 

Theorem 4.3. Let K be non-trapping for P and assume that there exits a pseudo-convex neigh- 
borhood ft D K of K. Assume that T is injective on £'(K). Then there exist k > and e > so 
that for any w\, W2, e-close to w\ and W2 in C k (ft), the corresponding operator I is still injective, 



and the estimate (4-20) holds with a constant C independent of the particular choice of w\, w^. 

Proof. By [9j, see also |16l Theorem 26.1.14], under the assumptions of the theorem one can con- 
struct a parametrix E so that 

EP = Id + R, 

where R has a smooth kernel. The parametrix E is not unique, even modulo smoothing operators. 
Loosely speaking, it is unique modulo smoothing operators if we fix an orientation on each connected 
set of bicharacteristics through K. The operator E has the mapping property E : HZ — > H s ~ l . 
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If we make P of order 1, then P would be microlocally equivalent to d/dx\; and then roughly 
speaking, E is integration w.r.t. xi in that representation in the direction of the chosen orientation. 



By (4.11 ), we have 



(4.23) 



EQ1 = Id + R, 



where Q is of order 1/2. Notice that E is of order 1, and X is of order —1/2. While the composition 



EQX a priori is of order 1 just based on the individual terms, it is actually of order as (4.23) 
shows. 

The construction of the Fourier Integral Operator (FIO) E is described in [16]. In order to get R 
above to be just of order —1, all microlocal constructions need to be done up to finite order only in 
order to satisfy finitely many symbol estimates, see, e.g., |15l Theorem 18.1.11'] and [26]. In each 
step, finitely many derivatives of the symbols are needed; therefore, finitely many derivatives of w\ 
and W2 are needed. Therefore, for some k, C k B (w\,W2) — >• R is continuous, where R : H s — > H s+1 
for a fixed s. 

The arguments below follow the proof of |29(, Proposition 5.1]. The idea is to correct the 
parametrix EQ by a finite rank operator so that the new Id + R would be injective. We should be 
able to do this because X is injective. 



Restrict equation (4.23) to K. In this stage of the proof, we will indicate the dependence on 
:= (wi,wi) by a subscript w. We can always assume that R w is self-adjoint because we can 
apply Id + i?^ to both sides of (4.23). The operator Id + R w has at most a finite-dimensional kernel 
V on L 2 (K). Since I w is injective on L 2 (K), I w : V — > T W V is an isomorphism; let B w be its 
inverse. Let also LT^ be the orthogonal projection to X W V . For w close to w as in the theorem, set 
B% :=E a Qu 1 + B V] Il w . Then 



(4.24) 



Bll^g = (ld + Rl)g, 



where RL := Rw + B w Il w X w is compact. We claim that Id + RL is injective for w 



w. Indeed, 

assume (Id + Rw)g = 0. Then (Id + R w )g + B w Il w I w g = 0. The first term is in V^; the second 
one — in V, therefore they are both zero. Thus g £ V, and B w Il w Z w g = 0. By the definition of 
B w and 14^, this implies g = 0. Therefore, Id + F$ w is injective, and actually invertible in L 2 (K). 
This property is preserved under small C k perturbations of w, k S> 1, as discussed above, with a 
uniformly bounded norm. The statement of the theorem now follows directly from (4.24). □ 



Theorem 4.3 and Theorem 4.2 imply the following generic uniqueness result. 



Corollary 4.2. Let K and O be as in Theorem \4-3\ For some k S> 1, there is an open dense set 
of pairs (wi,w<z) in C k (Q) so that the corresponding operator I is injective on £'(Vt), and satisfies 



the stability estimate (4-20) with a locally uniform constant. 



5. The non-linear Identification Problem 

Let (a, /) and (a, /) be two attenuation-source pairs. We will denote functions and oper ators 
related to (a, /) by placing a tilde over them. The difference v := u — u of the solutions of (2.1) 
solves 



(5.1) 

where 

(5.2) 



(9 ■ d x + a)v = 6f - u5a, v\e. x< i = B, 



5a :- 



Sf:=f~ f- 
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Therefore, 

(5.3) X a f - X a f = I w 5a + X a Sf, 
where I w is the weighted X-ray transform with weight 

(5.4) w = -e~ B ~ a u. 



of Proposition 3.1, as we should. 



We used here the obvious generalization of (2.2) for sources / dependent on 6 as well, see the 



remark following (2.3). If we replace a on the right with a, then we get the linearization formula 



5.1. A summary of the properties of the linearization 8X a j. We are in the situation of the 
previous section with 



(5.5) 



9i = $a, g 2 = 5f, wi 



W 2 



-Bd 



If 5a, 5f are given by ( |5.2[ ), then (4.7) is a non-linear equation, of course. If we treat them as 
independent (of a, a, f, f) functions then we have the linear problem that we analyzed above. 
Then 

-Bd-BJd-i 



(5.6) 



W 



l W , Wo := 



Ju), 



see (4.5). The characteristic variety £ in this case is given by 

(5.7) ^ = {u(x,em) = u(x,-e/\t\)}. 

The Hamiltonian po is then given by ( |4.12| ). Since an elliptic factor does not change the zero 
bicharacteritics, just their parameterization, the zero bicharacteristics are then given by the follow- 
ing Hamiltonian 

(5.8) H(x,0 = (u(x, 9) - u(x, -6)) \ e= ^ m - 
Recall that u us the solution of (|2.1l). 









* 



supp/ 



Figure 1. The zeros (x,8) of W are characterized by the property that the atten- 
uated integrals of / from x in the directions 6 and —0 are equal. The conormals £ 
to such 9 are the characteristic ones. If d g ±Wo(x, 6) / 0, then there is a smooth ray 
through x tangent to 8. 



We will summarize the properties of the rays, see Definition 4.1 in this case. Let X be the linear 
operator defined in (4.1) with weights w\ and W2 as in (5.5) but g\ = 5f and gi = 5a considered 
as independent functions. Notice that the rays depend on a and / only. On the other hand, a, f 
affect the weights in X. 
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For any x there is at least one ray through it which might be a point. 

The rays may not be smooth. Given (x,9) G R n x S , there is a smooth ray trough x in 

the direction of 9 if and only if Wq(x, 9) = and d g ±Wo(x, 9) / 0. 

Since 9 ■ d x Wo = 2f, the condition f(x) ^ is sufficient for P to be of real principal type 
at (x,g). 

A compact set K C R 2 is called non-trapping, if all rays eventually leave K. 

If K is non-trapping, then X, restricted to K, has a finite dimensional kernel, smooth enough 

if Ba and u are smooth enough near K. Also, (4.19) holds. 



• If X is injective on K, then it is stable, as well, with a loss of one derivative, i.e., (4.20) 
holds. If in addition K has a pseudo-convex neighborhood, then the injectivity is preserved 
under a small enough perturbation with a uniform stability estimate (4.20). 

• If Wq(xq, 9) = implies dg±Wo(xo, 9) ^ for all 9, then X is injective (and stable) restricted 
to functions supported in some neighborhood of xq- 

• If K is non-trapping, and Ba and u are analytic in a neighborhood of K, then X is injective 
(and stable). 

Remark 5.1. One important improvement in this case is to use the ellipticity of the second term X a 
in 5X a j, see Proposition 3.1 Let supp 8a C K%, supp 5f C K 2 , with K\ t % compact sets. Then 5X a j 
is elliptic on the set K~2 \ K\ because there (i.e., for 5a, 5f supported there), 5X a j(5a, 5f) = X a 5f. 
Therefore, for re cove ry of singularities we only need K\ to be non-trapping. If, in addition, 
on K\ \ K2, see (2.1 ), then it is enough to ask K\ n K2 to be non-trapping. 

5.2. Uniqueness and stability results. Our first main result about the identification problem 
is the following theorem. Recall that the requirement on Q to be pseudo-convex implies that K is 
non-trapping. 

Theorem 5.1. Let K C R 2 be a compact set and let Q D K be open. Let do, fo be of compact 
support so that their beam transforms Bao and Bfo are in C k (ft x S* 1 ). Assume that f2 is pseu- 
doconvex w.r.t. the Hamiltonian H defined in (5.8), related to ao and /q. Let do, fo be such that 
5X ao j , see Proposition 3.1, is injective on K. 

ao G C k and fj 



{a%, fx), (a 2 ,/2) with a,- 
(5.9) 



Then if k S> 1, there exists e > so that for any 
fo G C k supported in K satisfying 



\B(aj - a )\\ ck ^ xS i- ) + \\B(fj - /o) Hc^nxS 1 ) 



<£, j 



1,2, 



there exist constants C > 0, /i G (0, 1) so that 

(5-10) \\ax - a 2 \\ L 2 {K) + ||/i - f 2 \\ L ^K) < C\\X a Jx 



X ai fx 



l H 1 /2(Z)' 



Proof. By (5.3), we have 
(5.11) 



X a2 f2 — X ai fx — —I e -Ba lui 5a + X ai Sf + R 



where 
(5.12) 

and 5a = a2 — ax, 5f = / 2 — fx- Next, 

\\R\\l°° < C\\B(a 2 - ax) \\L°°(K)\\ a 2 



5X aiJl (5a,5f) + R, 



R — I(e- Ba i-e- Ba 2)m $a + X a2 ^ ai 5f, 



2)u% 



(5.13) 



a l\\L°°(K) + C\\B(a 2 - a>l)\\L°°(K)\\f2 - fl\\L°°(K) 



< C [\\5a\\ 2 LX{K) + \\5ff LOO {K) 



where C depends on an a priori bound of ||/i||L oo (n) which can always be found depending on ao, 
fo, e; by (p0). 
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We will apply [271 Theorem 2]. Set A(a, /) = X a f. Set also B x = L°° (K) x L°° (K) , B 2 = L°° (Z) . 
Then A : B\ — > B 2 is continuous. By (5.11 ) and (5.13), A is differentiable at (ai, /i) with a quadratic 
estimate of the remainder. 

By assumption, -BaolnxS 1 an d -^/olnxs 1 are m The same holds for the solution uq of (2.1) 
related to ao, /o- Since aj — oo £ Cq(K) and /j — /o G C$(K), we also get the same for i?Oj, 

1, 2 are O(e) perturbations of -Bao, 
4.3 to conclude that JX^jj is still 



J: 



1,2. Moreover, by (5.9), Baj, Bfj, and Uj, j 



and u 

Bfo, and «o in C k (0, x For k 3> 1, we apply Theorem 
injective, satisfying a stability estimate (4.20) with a constant C independent of ai, /i. Take s > 1 

3/2, to get 



in (4.20), for example, s 



1 5a I 



L°°(K) 



+ \\Sf\\ 



L°°(K) 



< 



C\\SX aiJl (5a,5f)\\ H s {z) . 



Based on that, we set 

B[ = B x = L°°{K), B' 2 = H 3 (Z). 
Then we have the following interpolation estimate 



<C\\h\\ 



i 

2 I 
L-'l 



1 

2 



where B2 



H (Z). We have all conditions met to apply |27l Theorem 2]. We therefore get that of 
k ^> 6 (k needs to satisfy both k > 6 and the requirements of Theorem 4.3), the stability estimate 
KM holds with ii = l/2. 



□ 



Remark 5.2. It is enough to assume that aj and fj, j = 0, 1,2, satisfy the regularity assumptions, 
instead of Baj, Bfj but that would be more restrictive. 



Remark 5.3. The support conditions for aj 
Remark 15.11 



ao and — fo can be relaxed to some extent as in 



Remark 5.4. The value for /j, that we got is fi = 1/2 but that was based on specific, and a bit 
arbitrary choice of the interpolation space H 6 . As shown in [271 Theorem 2] , and as can be easily 



seen from the proof, we can choose any [i > 1 in (5.10), as close to 1 as we wish, at the expense of 
increasing k. 

Next corollary is a generic local uniqueness and stability result on non-trapping sets. 

Corollary 5.1. Let K C R 2 be a compact set and £2 D K be open. Let ao, fo be of compact 
support so that their beam transforms Ba$ and Bfo are analytic in £1 x S 1 . Assume that Q is 
pseudoconvex w.r.t. the Hamiltonian H defined in (5.8), related to ao and /q. Then the conclusions 



of Theorem 5.1 hold 



This corollary implies in a trivial way also local uniqueness, in K, near a gene ric (dense and 
open in C k , k 3> 1) set of (a, /). The proof follows immediately from Theorem 



4.2 



The second corollary below states local uniqueness and stability in a small enough non-trapping 



set. 



Corollary 5.2. Let xo G R 2 is such that Ba, Bf are smooth near xq, and Wq satisfies (4-22) 
Then there exists an open set U B xq, so that for any K C U the conclusions of Theorem 5.1 hold 



Proof. By Corollary 4.1 if U is small enough, SX a j is injective on any compact set K C U. Then 
we apply Theorem 5.1 □ 
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5.3. Conditions for smoothness and analyticity of Ba and Bf. The results above require 
Ba and Bf to be either smooth enough or analytic in some open set fi. The smoothness, for 
example, certainly hold if a and / are smooth enough in Q but this is too restrictive. The following 
condition is sufficient. 

Proposition 5.1. Let Q, C R 2 be open. Let {cj}^ =l be a finite number of C k (respectively analytic) 
non-intersecting curves in R 2 \ Q so that a and f are C k /analytic in R 2 \ {cj}, up to the boundary 
on either side of each curve. Assume that each line through intersects every Cj transversely. 
Then Ba, Bf are in C k (£l x S 1 ), respectively analytic in $7 x S 1 . 

Proof. Near each (xq,9q) el!x S , Ba, and similarly Bf, is given by 

N pa +i(x,8) 

Ba(x,9) = ^2j a(x + te,d)dt, 

where «o = 0, aN+i = °°> an d the rest of the ay's are determined by the intersection points of the 
ray x + t9 with the curves Cj. The statement now follows directly from this representation. □ 

Remark 5.5. We presented the condition above in a form suitable for applications. For C°°, respec- 
tively, analytic regularity of a, / in Q, x S , it is necessary and sufficient to assume that a and / 
have the same regularity in Q; and a, /, have no C°°, respectively analytic singularities, conormal 
to some line through Q. The necessity follows from standard properties of the Radon transform 
to recover conormal smooth or analytic singularities. This condition is sufficient, because of the 
standard relation between the smooth/analytic wave front set of Ba or Bf on one side; and the 
Schwartz kernel of B and a or /, on the other. We sill skip the details. 



6. Further microlocal properties of 1 
Take the Fourier transform of L w (p0- L , 9) w.r.t. p to get 

(6.1) j e~ iX PL w f(p6 ± , ±6) dp = j e- iXe± -y W (y, ±9)f(y) dy. 
Set f = Xe^, A > 0, to get 

(6.2) J e-W\(I w f)( p e±,±Z- L /\Z\)dp= J e-'y<w{y,±em)f(y)dy. 
Take the inverse Fourier transform of both sides to get 

(6.3) w*(x,±D ± /\D\)f = (2tt)- 2 / e^^^fm/m^/l^dP^, 

where w*(x, ^D- 1 - /\D\) is the ^DO with amplitude a(x,y,^) = w(y, ±1;-*- /\t;\). The principal 
symbol is w(x, i^- 1 /^!). If w = 1, one can see that we get C\D\I[Iif on the right, and / on the 
left, which is just one of the inversion formulas for L\. 
Apply the described operation to the equation 

(6.4) I Wl gi + 1^5-2 = 0, 



compare with (4.7). We get 

(6.5) wl(x, ±D ± /\D\)g 1 + w* 2 (x, ±D ± /\D\)g 2 = 0. 



This is actually a system, see also (1.3). 



16 



PLAMEN STEFANOV 



Proposition 6.1. Let wip be two smooth weight functions, and let g = (51,(72) G £'(R n ). Then 
(6.6) Ig G H S {Z) 

if and only if 



(6.7) 



w* 1 (x,D ± /\D\) w* 1 (x,D ± /\D\) 
w* 2 (x,-D ± /\D\) w* 2 {x-D^/\D\) 



g G H S ^ 2 (R 2 ). 



Proof. Assume that the l.h.s. of (6.7) is in H s 1 > 2 . Then the r.h.s. of (6.3) with I w g replaced by 
Ig := I Wl gi + I W2 92 belongs to the same space. Take the Fourier transform of that to get, see also 
d6lb, 



(6i 



(A) 



s-1/2 



/ e-' lXp lg(p9 ± , ±0) dp G I? (R+ x S 1 , A dA dO) . 



Since the relation above holds with either choice of the ± sign, we can fix the positive one, and 
allow A to be negative, as well. Therefore, (\} s ~ 1 / 2 \\\ 1 / 2 J 7 ph - > .\I(p6 ± ,6) G L 2 (R x S 1 , dXd9). This 
easily implies, see e.g., the proof of [20l Theorem II. 5.1], that (X} s T ph ^xl(p9 ± ,9) G L 2 (R x S 1 ), 
which yields (6.6). 

Now, assume (6.6). Rev ersing the arguments above, we get (6.8). Take inverse Fourier transform 
w.r.t. £ = X6 1 - to get ((ojl). □ 



Proposition 6.1 reduces the problem of the microlocal invertibility of the FIO X to that of the 



matrix valued ^DO in (6.7) with a principal symbol 
(6.9) 



w;i(:r,#) W2(x,8) 
wi(x,—6) w 2 (x,-8) 



The determinant of the latter is W(x,£ /|£|), see (4.4) and (4.5). An immediate consequence of 
(6.7) is the following. For s ome matrix val ued classical ^DO P with a scalar principal symbol 
Po(x,^) = W(x, C^/ICI), see (4.12), relation (6.6) implies 

(6.10) 



This also follows from Proposition |4.1 



Pg G # S - 1/2 (R 2 ). 



Assume now that (4.18) is satisfied for some (xq, Oq). Then we can solve the equation W(x, 9) = 
for 9 G S 1 locally to get a smooth function 9{x). Since W is an odd function of 9, the same thing 
applies near the point (xo, — 9q), as well, with a solution —9{x). This implies that in a conic 
neighborhood of (xo,±9q) G S, the characteristic manifold £ is given by ^ _L /|^- L | = ±9(x). Set 
v ±(x) = Wj ( x , ±9(x)), j = 1, 2. Then Wj (x, ±£/\£ ± ) 



- v- 



(x) vanishes on S, and is therefore locally 
given by po(x,£) times a smooth function, homogeneous of order in £, hence a symbol. This 
implies that (6.7) can be written as 

(6.11) 



vt(x) Vo(x) 



g + (Q Q P + Q-i)g£ H^ixo,^ 



where Qq and Q—\ are classical ^DOs of order and —1, respectively. Using (6.10), we get 
(6.12) vfgi + vfg 2 + Q%g G H s ^ 2 (x ,^), 

with Q_ x of order —1, and the equations with the + and the — sign are actually linearly dependent 
up to the lower order term (including the possibility that one of them has zero coefficients). Now, 
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if the assumptions of Theorem 
we get 

(6.13) 



4.1 



are satisfied, (6.6) yields g G H s 3 / 2 . Then Q_ x g G H s 1 / 2 , and 



vfgi + vfgz G H'-V\xo,€q). 



Since the matrix in (6.11) has rank 1, only one of the equations (6.13) is relevant. This is an 



improvement over the estimate (4.20), that asserts that (6.6) implies g\ 2 G H s 3 / 2 , if suppg is 



supported in a non-trapping compact set. This improvement applies to the linear combination 



(6.13) only. 



6.1. Applications to the linearized Identification Problem. Let X 

tion of X a f, see Proposition 



3.1 



Then /,-, Wj are given by (5.5). 



8X a j be the lineariza- 
The determinant W can be 



replaced in the analysis by Wo, see (5.6). Notice that W2 > 0. The discussion above yields the 
following. 

Proposition 6.2. Fix (x o ,0 o ) G R 2 x S 1 . Let W (x ,6>o) = and d 9 ± W (x , O ) / 0. Let 
(6.14) v(x) = u(x, 0(x)), for x near xq, 

where 9{x) is the unique local solution of Wq(x,9) = with 9(xq,9q) = 9q, and u is defined by 
(2.2). Then, if the assumptions of Theorem 4-1 are satisfied, and if 8X a j(8a,8f) G H s , we have 



(6.15) 



,8a- 8 f G H s -^ 2 (x ,±e^). 



Proof. In this particular case, w\ = —e Ba u, W2 = e Ba . Under the non-degeneracy assumption 
on Wo, W2 > 0, and w\ = —UW2- Divide by the elliptic factor u>2 in either of the two relations 
(67131) to get KTEh. □ 



Remark 6.1. Theorem 4.1 says that under the non-trapping condition we can recover WF{j ^2) 
with a loss of one derivative, compared to the standard X-ray transform. On the other hand, 



Proposition 6.2 says that under the additional mild condition on W, one can recover the wave 



front set of the linear combination (6.15) without loss. This has the following implications for the 
recovery of a and /: we can expect v 8a — 8f to be recoverable in a more stable way than either 8a 
or 8f. 



Remark 6.2. We need to assume that the assumptions of Theorem 4.1 ar e satis fied just to conclude 
that / G H s - 3 / 2 ; and then to deduce that Qt x f G H 8 - 1 / 2 , see (6.12) and (6.13). If we know a priori 
that / has certain regularity, then we can use that fact instead. In applications, it would be natural 
to assume that (8a, 8f) G L? . Let us assume that the measurements show that 8X a j(8a, 8f) G H 3 / 2 
(or better). Then we conclude that v8a — 8f G H 1 , that in particular excludes jump types of 
singularities at smooth surfaces of that particular linear combination. There is no need to assume 
the trapping condition for this argument. 



7. The radial case 

As explained in the Introduction, the thorough study of the case of radial a and / is behind 
the scope of this work. The purpose of this section is to present a case, where the rays can be 
easily computed, when both the linearized map, and the non-linear one have huge kernels if the 
non-trapping assumption is not satisfied. So at least in the cases described below, the non-trapping 
assumptions is not only sufficient but also necessary for the problem to be "well-behaved" . 
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7.1. The linearized map for a simple radially symmetric example. We start with perhaps 
the simplest example. Let 1b(o,i) be the characteristic function of the unit disk. We study the 
linearization 5X w.r.t. (a, f) near 

(7.1) a = 0, / = 1b (0 ,i). 



We will choose perturbations of those a and / supported in B(0, 1) only. The weight w, see (3.2) 
or (3.3), restricted to the unit disk, is given by 



(7.2) 



w(x,e) = -Ji-(e ± -x) 2 



■ X. 



Then, see (5.6), 
(7.3) 



-26 -x. 



The Hamiltonian H, up to a constant factor, is as in Example 4.1 Indeed, by (5.8), H 
f-L/lfl = -2(xi6 - z 2 fi)/|f|. Therefore, -2\£\H is the symbol of 

xiD 2 - x 2 Di = -id/d(j>, 



-2x 



where 4> is the polar angle in the x space. The bicharacteristics are given by (4.16). In particular 



the rays are the concentric circles |x| = R, R > 0, including the degenerate case x = 0. As before, 
K C B(0, 1) is non-trapping, if and only if K does not contain an entire circle of that kind, see 
Figure [2j 





Figure 2. The rays of Example 4.1 in the unit disk and an example of a non- 
trapping K, left; and a trapping K, right. 

The equation 6X(5a, Sf) = can then be written as 

-/ (jl-(6 ± -x) 2 + 9-x)6ads+ [ 6fds = 0, 

where l z g is the line through z G 9 1 - in the direction of 6, and ds is the natural measure on it. The 
integral over the line I z -q would produce the same term with 9 ■ x replaced by —9 ■ x. Therefore, 
both the even and the odd part w.r.t. 9 above vanish: 



(7.4) 



f Jl- (9 ± -x) 2 5ads+ [ 5fds = 0, [ 9-x5ads = 0. 



The third integral is the X-ray transform of the vector field {5a)x. It is well known that we can 
only determine the curl of that, i.e., 

(x\d2 — X2d\)5a = 0. 
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In other words, 5a needs to be radial. Then the first term in (7.4) is invariant under rotations of 
(x,6), i.e., when we consider (x,6) as points in the unit tangent bundle. Then so is the second 
term. Apply /[ to it, and we get that \D\~ 1 5f is radial, as well. Then so is Sf. 

Therefore, the kernel of 5X(5a,5f) consists of radial 5a and 5f that are connected by the first 
identity in (7.4). Since the weight there is constant along the lines, using Radon transform notation, 
Rh(p, uj), we get 



(7.5) 



y/l - p 2 R5a - RSf = 0. 



It follows from the analysis below that there exists an infinite dimensional space of pairs (5a, 5f) 
satisfying (7.5). Indeed, for any radial 5a £ C^°(B(0, 1)), we can solve ( |7.5[ ) for 5f, and vice versa. 

Going back to ( 7.4 ) , th e arguments in the proof of Proposition |6.2| (or its conclusion directly, 
together with Remark 6.2) show that 5X(5a,5f) 6 H s and (5a, 5 f) € _£P -3 / 2 imply 



.r 



l 5a -5f eH s 



-1/2 



in the annulus U := {0 < |x| < 1}; i.e., the singularities of that particular linear combination in U 
can be recovered without a derivative loss. Note that for any x £ U, the characteristic directions 
(zeros of W) are given by 8 = ±x /\x\, and the characteristic codirections — by £ = ±x/|x|. Then 
the integral of /, starting from x, in a characteristic direction 9 is exactly \J\ — |x| 2 |. This is the 
value of u for characteristic directions, see (5.8) and (6.14), and confirms (6.15). 



7.2. The linearized map for a = and / radial has an infinite dimensional kernel. Let 

now a and / be general radial smooth functions of compact support. Then the characteristic variety 
of Example 41 and Section [7] 



So = {(:#,£); x and £ are collinear} 

is included in the characteristic variety £ in this case but the latter can be larger. The Hamiltonian 
curves (with rays |x| = R > 0) of those examples are still Hamiltonian curves in the present case 
but there may be more. If / > in -B(0, 1), then it is easy to see that the Hamiltonian flow over 
B(0, 1) is the same. 

We study now 5X a j for 

(7.6) a = 0, / radial. 

We also assume that / is smooth and has compact support. With some abuse of notation, we 



replace / by / = /(|x|), where / has even smooth extension. By Proposition 3.1 
(7.7) 5X 0J (5a, Sf) = -I JBf Sa + I 5f, 



see (4.6). We restrict SXqj to radial 5a, Sf, as well. 



We will use Radon type of parameterization for IjgjSa by setting oj = . Write 
RjBfSa(p,uj) = IjBf5a(puj,—uj i ~)= I 5(p — uj ■ x)Bf(x,uj ± )5a(x) dx. 



Here 5 is the Dirac Delta function, not to be confused with the variation symbol in 5a, Sf. Since 
/ is radial, for any rotation U, we have Bf(Ux, Uuj- 1 ) = Bf(x,oj^~). Since 5a is radial as well, we 
easily get that IjBfSa is independent of oj, i.e., IjBfSa = IjBf5a(p). We claim that IjBf5a(p) is 
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an even function of p. Indeed, set oj = (1,0). Then 

RjBfH~P) = J 5 (~P - xi)Bf(x, (0, l))Sa(x) dx 

= J 5(p + x\)Bf(x, (0, l))5a(x) dx because 5 is even 

= J 8(p — x\)Bf(x, (0, l))8a(x) dx after the change x\ \-¥ —x\ 
= R JB f8a(p). 

In the last equation, we also used the fact that / is radial. 
To study the kernel of 5X j, we write, see (7.7), 

(7.8) R JBfS a-R8f = 0, 

where, with some change of notation again, R is the classical Radon transform acting on radial 
functions, i.e., considered as a map on functions of a single variable. It is easy to see that, see also 
H2I, 



Rg(p) = 2 J g(y/p 2 + t 2 )dt, p>0. 



It is known, see [20J, and can be easily seen that this equation can be written in the form 

Rg(p) = 2^ (l-^T) 1 9(r)dr, p>0. 

This an equation of Abel type with explicit inversion given by (see [lOl [20] ) 

1 f'°° c\ 

(7.9) g(r) = — {p 2 -r 2 )- l / 2 -Rg(p)dp. 

TT J r dp 

Moreover, the Abel transform R is given by a composition of the cosine Fourier transform F c and the 
zero order Hankel one Hq (see |10j). with a proper normalization, i.e., R = F c Hq. If h € C°°(R+) 
is of compact support, and admits a smooth even extension, then we get a direct confirmation that 
the equation Rg = h has a (unique) solution given by g = HoF c h. Indeed, for such h, F c h has 
smooth even extension in the Schwartz class, and then HoF c h is well defined and solves IZg = h. 
This shows that the function 5 in (7.9) is given by 

(7.10) S f(r) = -~J (p 2 -r 2 )- 1 / 2 -Ij Bf 5a(pu J ,u ± )dp, 

see also (7.7). We recall that IjBfSa is independent of oj. We summarize this into the following. 

Proposition 7.1. Let f G C^°(R 2 ) be radial. Then the linearized map 5Xqj (with a = 0) has an 
infinite dimensional kernel, including all radial pairs (5a, 5f) with 5a smooth function of compact 
support, and 5f given by (7.l0\j . 

In other words, besides the inability to recover the singularities (without support restrictions), 
we actually have an infinite dimensional kernel. Therefore, in this case, the non-trapping condition 
is a necessary condition for the problem to be well posed, as well. 
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7.3. Non-uniqueness for the Identification Problem for radial a, / near a = 0. We show 
next that not only does the linearized map SX a t can have an infinite dimensional kernel in the 
case above, but the non-linear map (a, /) i— > X a f has a rich set of radial pairs with the same image. 

Theorem 7.1. Let a G Cq° and f G Cq° be radial. Then there exists fo G Cq° so that 

(7.11) X a f = X fo. 

Proof. We will work again with the Radon transform parameterization R a f(p,ui) = X a (puj, —uj 1 ) 
instead, see (2.5). As above, it is straightforward to check that 

Rafip, w) = R a f(-P, -w). 

We saw above that Rof(p, uj), denoted there by Rf, is actually independent of uj, and an even 
function of p. Then for any k = 0,1, . . . , 



J R a f(p, u)p h dp = C k = const., 



and Cfc = if k is odd. Therefore, the integral above is a restriction of the homogeneous polynomial 
Cfc|^| fc to the unit sphe re. Therefore, R a f G Sh, and by the Helgason range characterization 
theorem, see [12J, (7.11) holds with some /o € 5(R 2 ). By the support theorem, /o is compactly 
supported. □ 



We can actually make this constructive. By (7.9), writing /o = foi r ), we get 

1 f°° H 

fo(r) = — ip 2 ~ r 2 )- l ' 2 -X a fipu, -uj 1 ) dp, 

7T J r dp 

recall that X a f{puj, ui 1 ) is independent of uj. 

Appendix A. I£I a as a ^DO 

As explained in Section [2j we view X a f and I w f as functions on Z, with a natural measure dz 
there. Then X a , and more generally, I w have well defined transpose (w.r.t. the distribution pairing) 
and conjugate (w.r.t. the L 2 product) operators X' a and X*; and I' w , 1^, respectively. Below, we 
use the notation 8± for the line given by s i— y pO 1 - . 

Proposition A.l. 

I t w i/>(x)= I w(x,6)ij(x- (x- 6)9,9) d9. 



Proof. For G C§°(R 2 ), ijj G Cg°(Z), we have 



(I w (j))^dz= / / wiz + s9,9)(piz + s9)^iz,d)dsdzd9. 

J Z JR 

Set x = z + s9, z G 9±. For any fixed 9, (z, s) t— )■ x is a diffeomorphism with a Jacobian equal to 1. 
Its inverse is given by 

z = x — (x ■ 9)9, s = x ■ 9. 

Therefore, 



L 



(I w <t>)ip dz = / / w{x,0)<f>{x)ij}{x- (x- 9)9,9) dxd9, 
and this proves the proposition. □ 
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Proposition A. 2. For any two smooth functions a and b, 

l' b I a f(x) = / _' y| /(y)dy, 

w/iere 

(A.l) y, 0) = a(x, 0)6(y, 0) + a(x, -0)b(y, -9). 

Moreover, I' b I a is a classical ^DO of order — 1 with amplitude 

(a.2) ^ (afo^/KD&foeVif i) + o(x, -e^/iei) -e7lf I)) 

and principal symbol 

^( a (x,e/\mxA ± M)+a(x-i; ± /\mx-!i ± m) 



Proof. By Proposition ! A.l 



l' a I b f(x) = J b[x, 9) j a(x - (x ■ 9)9 + W, 9)f(x - (x ■ 9)9 + td) dt d9 
b(x,9) [ a(x + t9,9)f(x + t9)dtd9. 



Split the i-integral in two parts: for t > and for t < 0, and replace t by — i in the second one to 
get 



l' a I b f( x ) = J b(x,9) J a(x + t9,9)f(x + t9)dtd9 

P f'OO 

(A.3) = / b(x,9) / a(x + t9,9)f(x + t9)dtd9 

Js 1 Jo 



s 1 Jo 

P f'OO 

+ / b(x,9) / a(x-t9,9)f(x-t9)dtd9. 
Js 1 Jo 

Replace —9 by 9 in the second integral to get 

(A A) l' a I b f(x) = / [b(x,9)a(x + t9,9) + b(x,-9)a(x + t9,-9)]f(x + t9)dtd9. 

Js 1 Jo 

Pass to polar coordinates x = y + t9, centered at y to finish the proof. 

To write as a ^DO, recall that if the Schwartz kernel of a linear operator is given by 
K(x, y, (x — y)/\x — y\), then it is a formal ^DO with an amplitude given by the Fourier transform 
of K w.r.t. the third variable. Therefore, I*Ib is a formal ^DO with amplitude 



e 12 " 



\z\- l A(x,y,z/\z\)dz= [ e irH A(x, y, 9) dr d9 = vr / A(x, y, 9)5(9 ■ d9 

Jn+xs 1 Js 1 

7T 

= W\ 

We used here the fact that A is an even function of 9 and that the inverse Fourier transform of 1 
is S, see also |14l Theorem 7.1.24]. Since this is a homogeneous function of £, with an integrable 
singularity that can be cut-off resulting in a smoothing operator, this completes the proof. □ 



A(x,y,em\) + A(x,y,-^m). 
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The mapping properties of those operators are well understood even in the more general setting 
of the weighted geodesic transform. We summarize those properties below. Recall the definition 
of the Sobolev space H S (Z) in (2.6) first. Given a compact set K C R 2 , we also use the notation 
H S (K) to denote the closed subspace of the distributions in H S (R?) supported in K, see [32l 
Chapter 4.5], where those spaces are denoted by H S K (M). 

Proposition A. 3. For any compact set K C R 2 , and any s > 0, 

(A.5) I w :H s -V\K)^H s comp (Z), l' w : H^Z) -> Hf oc (B?) 

are continuous. 

Proof. We follow the proof of Proposition 5.1 in [30 1. We can always assume that w is extended 
smoothly for x outside K so that it vanishes outside a small neighborhood of K. Then we can 
replace R 2 and Z by compact manifolds, as explained in Section [5J and work with / £ C°°(T 2 ). 
Note first that 1'^ : H s — > H s+1 . Next, if s > is an integer, for / supported in K, 



(A.6) \\I w ff H s(z) < C E | {9iI w f,I w f) L2(z) =C7£ (l^ p I w f,f) 

j<2s j<2s 



The term I^dpl w f is a sum of weighted X-ray transforms of derivatives of / up to order 2s + 1, 
and is therefore a \I/DO of order 2s. This easily implies that for / £ C°°(T 2 ), 

\\Iwf\\ 2 H s( Z ) < C\\f\\ HS _ 1/2 ^ K y 

The case of general s > follows by interpolation. The estimate then holds for any / £ fP _1 / 2 (T 2 ), 
and therefore, for any / £ H s ~ l l 2 {K), as well. 

To prove the second estimate, notice first that d a I^ip is a sum of operators of the kind 1^ but 
with different weights applied to p-derivatives of tp up to order \a\. Then for any integer j > 0, 

lifXWv] = \(U,d^) L2{z) \ < C\\f\\ L2 U\\ HJ - 1/2 . 
This proves the second estimate for s = 0, 1, . . . . For general s > we use interpolation. □ 
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